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ABSTRACT 



in 



We use a combination of X-shooter spectroscopy, ULTRACAM high-speed pho- 
tometry and SOFI near-infrared photometry to measure the masses and radii of 
both components of the eclipsing post common envelope binaries SDSS J121258.25- 
012310.1 and GKVir. For both systems we measure the gravitational redshift of the 
white dwarf and combine it with light curve model fits to determine the inclinations, 
masses and radii. For SDSS J1212-0123 we find an inclination of i = 85.7° ± 0.5°, 
masses of M WD = 0.439 ± O.OO2M and M scc = 0.273 ± O.OO2M and radii R wu = 
0.0168 ±0.0003R Q and R scc = 0.306 ±O.OO7R . For GKVir we find an inclination of 
i = 89.5° ± 0.6°, masses of M WD = 0.564 ± O.O14M and M scc = 0.116 ± O.OO3M 
and radii R WB = 0.0170 ± O.OOO4R and R scc = 0.155 ± O.OO3R . The mass and 
radius of the white dwarf in GKVir are consistent with evolutionary models for a 
50, 000K carbon-oxygen core white dwarf. Although the mass and radius of the white 
dwarf in SDSS J1212-0123 are consistent with carbon-oxygen core models, evolution- 
ary models imply that a white dwarf with such a low mass and in a short period 
binary must have a helium core. The mass and radius measurements are consistent 
■ with helium core models but only if the white dwarf has a very thin hydrogen envelope 

(Mh/Mwd ^ 10~ 6 ). Such a thin envelope has not been predicted by any evolutionary 
models. The mass and radius of the secondary star in GKVir are consistent with evo- 
lutionary models after correcting for the effects of irradiation by the white dwarf. The 
secondary star in SDSS J1212-0123 has a radius ~ 9 per cent larger than predicted. 

Key words: binaries: eclipsing - stars: fundamental parameters - stars: late-type - 
white dwarfs 



1 INTRODUCTION 

Detached eclipsing binaries are a primary source of accu- 
rate physical properties of stars and stellar remnants. A 
combination of modelling their light curves and measur- 
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ing the radial velocities of both components allows us to 
measure masses and radii to a precision of better than 
1 per cent (e.g. Andersen _ 199ll: Southworth et al.1 120051 : 



ISouthworth et al.ll2007l ; lTorres et al.ll20ld ). These measure- 
ments are crucial for testing theoretical mass-radius rela- 
tions, which are used in a wide range of astrophysical cir- 
cumstances such as inferring accurate masses and radii of 
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transiting exoplanets, calibrating stellar evolutionary mod- 
els and understanding the late evolutio n of mass transferring 
binaries such as cat aclysmic variables dLittlefair et al.|[2008l ; 
ISavourv et al.ll201 ll ). Additionally, the mass-radius relation 
for white dwarfs has played an impor tant role in estimat - 
ing the distance to globular clusters l|Renzini et al.l | l996l) 
and t he determination of the age of the galactic disk l|Woodl 
Il992h . 

On the one hand, although ubiquitous in the solar 
neighbourhood, the fundamental properties of low-mass M 
dwarfs ar e not as well und erstood as those of more mas- 
sive stars (jKraus et al.l201lh . There is disagreement between 
models and observations, consistently resulting in radii up to 
15 per cent larger and effective temperatures 400K or more 
below theoretical predictions ijRibad l20od : iLopez-Moralesl 

2007) . These inconsistencies are not only seen in M dwarf 

J II 7i 1 II 1 

eclipsing binaries (Baylcss & Orosz 2006: Kraus et al. 2011) 
but also in field stars (Be rger et al] 120061 : iMorales et al.l 

2008) and th e host stars of transiting extra-solar planets 
l|Torresl [20071 ). 

On the other hand, the mass-radius r elation for white 
dwarf s is all but untested observationally. IProvencal et al.l 
(1 19981 ) used Hipparcos parallaxes to determine the radii 
for white dwarfs in visual binaries, common proper-motion 
(CPM) systems and field white dwarfs. However, the ra- 
dius measurements for all of these systems still rely to 
some extent on model atmosphere calculations. For field 
white dwarfs t he m ass determinations are also indirect. 
iBarstow et~ai1 (|2005l ) used Hubble Space Telescope/STIS 
spectra to measure the mass of Sirius B to high precision, 
however, their radius constraint still relied on model at- 
mosphere calculations and is therefore less direct when it 
comes to testing white dwarf mass-radius relations. Dou- 
ble white dwarf eclipsing binaries potentially allow ex- 
tremely precise measurements of white dw arf masses and 
radii but have only recently been discovered ( Stcinfadt et al. 
20ld:IParsons et alJl^Tlal : iBrown et al.ll201ll ; IVennes etall 



2011bl ). To date, only a handful of white dwarfs have 
had their m asses and radii model -independe ntly measured, 
V471 Tau ijO'Brien et all l200ll). NNSer [IParsons et all 
l2010al ). SPSS J0857+0342 (iParsons et alJlioilbl ) and SPSS 
J1210+3347 jPvrzas et ai1l201ll ). All of these systems are 
short period post common envelope binaries (PCEBs), 
demonstrating the potential of these systems for testing 
mass-radius relations of both low mass stars and white 
dwarfs. 

In this paper we combine X-shooter spectroscopy, UL- 
TRACAM high-speed photometry and SOFI near-infrared 
photometry to determine precise system parameters for 
the eclipsing PCEB s SPSS J121258. 25-012310.1 (henc eforth 
SPSS J1 212-0123) jNebot Gomez-Moran et al] 120091 ) and 



GK Vir (|Green et al. 19781 ). We then compare our mass and 
radius measurements with theoretical mass-radius relations 
for white dwarfs and low-mass stars. 



Technology Telescope (NTT) at La Silla in April 2010. 
ULTRACAM is a high-speed, triple-beam CCP camera 
ijPhillon et al.ll2007n which can acquire simultaneous images 
in three different bands; for our observations we used the 
SPSS u' , g' and i' filters. A complete log of these obser- 
vations is given in Table [T] The data collected for GK Vir 
were combined with previous observations of the system 
taken with ULTRACAM (see IParsons et all l2010bl for de- 
tails of these observations). We windowed the CCP in order 
to achieve exposure times of ~3 seconds which we varied 
to account for the conditions. The dead time between expo- 
sures was ~ 25 ms. 

All of these data were reduced using the ULTRACAM 
pipeline software. Pebiassing, flatfielding and sky back- 
ground subtraction were performed in the standard way. 
The source flux was determined with aperture photome- 
try using a variable aperture, whereby the radius of the 
aperture is scaled according to the full width at half maxi- 
mum (FWHM). Variations in observing conditions were ac- 
counted for by determining the flux relative to a comparison 
star in the field of view. Apparent magnitudes and coor- 
dinates for each of the comparison stars used are given in 
Table[2] We flux calibrated our targets by determining atmo- 
spheric extinction coefficients in each of the bands in which 
we observed and calculated the absolute flux of our targets 
using observations of standard stars (from lSmith et al1l2002l ) 
taken in twilight. Using our absorption coefficients we ex- 
trapolated all fluxes to an airmass of 0. The systematic error 
introduced by our flux calibration is < 0.1 mag in all bands. 



2.2 SOFI J-band photometry 

We observed both GK V ir and SPSS J1212-0123 with SOFI 
jMoorwood et al.lll998l ) mounted at the NTT in April 2010 
and April 2011. The observations were made in fast pho- 
tometry mode equipped with a J-band filter. We windowed 
the detector to achieve a cycle time of ~ 10-15 seconds and 
offset the telescope every 10 minutes in order to improve 
sky subtraction. A summary of these observations is given 
in Table Q] 

The dark current removal (which also removes the bias) 
and flatfielding were performed in the standard way. Sky 
subtraction was achieved by using observations of the sky 
when the target had been offset. The average sky level was 
then added back so that we could determine the source flux 
and its uncertainty with standard aperture photometry, us- 
ing a variable aperture, within the ULTRACAM pipeline. 
A comparison star was used to account for variations in ob- 
serving conditions, details of these are given in Table(2] Flux 
calibration was done using the comparison s tar J-band mag- 
nitud e retrieved from the 2MASS catalogue ( Skrutskie et al. 
120061) . 



2 OBSERVATIONS AND THEIR REDUCTION 
2.1 ULTRACAM photometry 

GKVir and SDSS J1212-0123 were observed with ULTRA- 
CAM mounted as a visitor instrument on the 3.5m New 



2.3 X-shooter spectroscopy 

GKVir and SDSS J1212-0123 w ere both observed with X- 
shooter jD'Odorico et al]|2006l ) mounted at the VLT-UT2 
telescope on the night of the 4th of April 2010. Due to the 
long orbital periods of both systems (~ 8 hours) we targeted 
the quadrature phases, since these phases are the most sen- 
sitive to the radial velocity amplitude. Conditions through- 
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Table 1. Journal of observations. Exposure times for X-shooter observations are for UVB arm, VIS arm and NIR arm respectively. The 
primary eclipse occurs at phase 1, 2 etc. 



Date at 


Target 


Instrument 


Filter (s) 


Start 


Orbital 


Exposure 


Conditions 


start of run 








CUT) 


phase 


time (s) 


(Transparency, seeing) 


nm n 1 c\ a /ntr 

2010/04/05 


C< T\ C 1 C T "1 O "1 o nl no 

bDbb J 1212-0126 


X-shooter 




00:42 


n rr t n 

0.57-0.79 


onn ooo inn 

300,338,100 


Excellent, ~1 arcsec 


nm n 1 c\ a /ncr 

2010/04/05 


GK Vir 


X-shooter 




04:28 


n 1 o n n/" 

U.12— U.2o 


onn ooo inn 

300,338,100 


Excellent, ~1 arcsec 


nm n 1 n a /nrr 

2010/04/05 


C T\ O C T "1 O "1 o moo 

bDbb J 1212-U12o 


X-shooter 




05:57 


n oo n 

0.23-0.45 


onn ooo inn 

300,338,100 


Excellent, ~1 arcsec 


on 1 n / r\ a /nrr 

2010/04/05 


GK Vir 


X-shooter 




07:55 


n rr o n 

0.53-0.75 


onn ooo inn 

300,338,100 


Good, ~1.2 arcsec 


nm n 1 c\ a /01 
2U1U/U4/21 


GK Vir 


T TT T^TU A n A AT 


, 

u g i 


07:18 


n nn i nrr 

0.90-1.05 


3.0 


Good, ^1.2 arcsec 


nm n I c\ a /on 
2U1U/U4/22 


C< T\ C 1 C T1 O "1 o m oo 

bDbb J 1212-0126 


T TT A n A AT 


„ i~ia> 
u g i 


02:26 


n a n n cr rr 

0.40-0.55 


2.9 


Good, ~1.2 arcsec 


2U1U/U4/22 


GK Vir 


T TT T^TJ A n A AT 


„ l ~l a< 
u g i 


04:14 


0.45-0.55 


3.0 


Excellent, ^1 arcsec 


om n /n/i / o o 

2010/04/23 


( ■ t "~\ oo TTOio n 1 n o 

bDbb J1212-0123 


T TT TITl A A H T 

UL1RACAM 


! 1 -1 

u g i 


02:46 


n a n n f p" 

0.40-0.55 


2.9 


Excellent, <1 arcsec 


om r\ /n,i / o o 

2010/04/23 


GK Vir 


T TT TF> A f~1 A TV \ 

UL1HAOAM 


u g i 


04:45 


n -i n n rr rr 

0.40-0.55 


3.0 


Excellent, <1 arcsec 


2U1U/U4/24 


O T"\ O O T1 oi o m oo 

bDbb J 1212-0126 


T TT TF> A f~* A TV \ 

UL1HAOAM 


u g % 


22:57 


0.90—1.15 


2.9 


Good, ~1.2 arcsec 


riAl n 1 r\ A /OCT 

2010/04/25 


GK Vir 


T TT T-Y) A A Tv /T 

UL1KAOAM 


„ l ~* A' 

u g i 


06:29 


n a n n 

0.40-0.55 


3.0 


Variable, 1-2 arcsec 


i"iai r\ 1 f\ a /on 

2010/04/29 


GK Vir 


SOFI 


■1 


01:16 


0.40-0.55 


15.0 


Excellent, <1 arcsec 


o m n / n /i / in 

2010/04/30 


GK Vir 


SOFI 


J 


02:17 


n a rr n t;n 

0.45-0.60 


15.0 


Good, ~1.2 arcsec 


n m n / n /i / o rt 

2010/04/30 


O T"\ O C 1 T1 01 o m oo 

bDbb J 1212-0126 


SOFI 


J 


03:56 


0.40-0.55 


10.0 


Excellent, ~1 arcsec 


n m n / n /i /on 

2010/04/30 


GK Vir 


SOFI 


■1 


06:16 


n nn 1 nrr 

0.90-1.05 


15.0 


Excellent, ^1 arcsec 


nm n /ntr /m 

2010/05/01 


GK Vir 


SOFI 


■1 


02:42 


0.40-0.55 


15.0 


Good, ~T.2 arcsec 


nm 1 /n/i /no 

2011/1)4/1)2 


O T"\ O C 1 TIOIO mOO 

bDbb J 1212-0126 


SOFI 


J 


01:07 


c\ a o n cn 


10.0 


Excellent, <1 arcsec 


nni -i /n /] /no 
2U11 /U4/U2 


ijrv Vir 


DtJr 1 


j 
J 


04:05 


n oo n en 


15.0 


Excellent, <1 arcsec 


nni -i /n ,i /no 
2U11 /U4/U2 


ijrv Vir 


DtJr 1 


j 
J 


08:02 


U.oD— l.Uo 


15.0 


Excellent, <1 arcsec 


om 1 /n/i /no 

2011/04/03 


cnno ttoio moo 

bDbb J1212-U123 


SOFI 


■1 


01:04 


n ,tn n ci 

0.40-0.61 


10.0 


Excellent, ^1 arcsec 


2011/04/03 


GKVir 


SOFI 


J 


05:03 


0.40-0.61 


15.0 


Excellent, ~1 arcsec 


2011/04/04 


SDSS J1212-0123 


SOFI 


.1 


01:34 


0.44-0.62 


10.0 


Excellent, ~1 arcsec 


2011/04/04 


GKVir 


SOFI 


J 


05:33 


0.37-0.66 


15.0 


Excellent, ~1 arcsec 


2011/04/05 


SDSS J1212-0123 


SOFI 


J 


01:04 


0.35-0.57 


10.0 


Excellent, ~1 arcsec 


2011/04/05 


GKVir 


SOFI 


J 


05:56 


0.95-1.20 


15.0 


Excellent, ~1 arcsec 


2011/04/05 


SDSS J1212-0123 


SOFI 


J 


05:06 


0.85-1.06 


10.0 


Excellent, ~1 arcsec 


2011/04/05 


GKVir 


SOFI 


J 


06:53 


0.43-0.69 


15.0 


Excellent, ~1 arcsec 


2011/04/06 


SDSS J1212-0123 


SOFI 


J 


01:54 


0.43-0.60 


10.0 


Excellent, <1 arcsec 


2011/04/06 


GKVir 


SOFI 


J 


03:19 


0.90-1.17 


15.0 


Excellent, <1 arcsec 


2011/04/06 


SDSS J1212-0123 


SOFI 


.1 


05:59 


0.94-1.12 


10.0 


Excellent, <1 arcsec 


2011/04/06 


GKVir 


SOFI 


J 


07:32 


0.42-0.62 


15.0 


Excellent, <1 arcsec 



Table 2. Comparison star apparent magnitudes and coordinates. Magnitudes listed are for those bands for which the comparison star 
was used. 



Target 


Comp 


Comp 


u' 


9' 


i' 


J 


star 


RA 


Dec 










GKVir 


14:15:22.86 


+01:19:12.7 


14.5 


13.5 


13.1 




SDSS J1212-0123 


12:12:54.97 


-01:20:25.0 


13.8 


12.5 


11.9 




GK Vir 


14:15:29.31 


+01:17:37.9 








13.7 


SDSS J1212-0123 


12:13:00.51 


-01:23:23.3 








13.4 



out the night were excellent with seeing consistently below 
1 arcsec. Details of these observations are listed in Table [1] 
X-shooter is a medium resolution spectrograph consist- 
ing of 3 independent arms that give simultaneous spectra 
longward of the atmospheric cutoff (0.3 microns) in the UV 
(the "UVB" arm), optical (the "VIS" arm) and up to 2.5 
microns in the near- infrared (the "NIR" arm). We used slit 
widths of 0.8", 0.9" and 0.9" in X-shooter's three arms and 
binned by a factor of two in the dispersion direction, result- 
ing in a resolution of R ~ 7, 000. We used exposure times of 
300 seconds in the UVB arm, 338 seconds in the VIS arm 
and 3 x 100 seconds in the NIR arm. After each exposure 
we nodded along the slit to help the sky subtraction in the 
NIR arm. 

The reduction of the raw frames was conducted us- 



ing the standard pipeline release of the X-shooter Common 
Pipeline Library (CPL) recipes (version 1.3.7) within ES- 
ORex, the ESO Recipe Execution Tool, version 3.9.0. The 
standard recipes were used to optimally extract and wave- 
length calibrate each spectrum. For the NIR arm we com- 
bined frames taken at different nod positions to improve 
the sky subtraction, however, this does result in a reduc- 
tion in orbital phase resolution. The instrumental response 
was removed by observing the spectrophotometric standard 
star GJ 440 and div iding it by a flux table of the same star 
l]Hamuv et al.lll992T ) to produce the response function, this 
was also used to apply a telluric correction to the spectra. 
We then heliocentrically corrected the wavelength scales of 
each of the spectra. We achieved a signal-to-noise (S/N) for 
GKVir of ~ 20 in the UVB arm per exposure, ~ 15 in the 
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Wavelength (microns) 

Figure 1. Averaged X-shooter spectrum of SDSS J1212-0123. The SDSS u'g'r'i'z' and 2MASS JHK magnitudes and filter widths are 
also shown. A zoom in on the white dwarf features are shown inset with the narrow Call and Mgn absorption features labelled, although 
the Can feature is somewhat filled in by emission from the secondary star. 



VIS arm per exposure and ~ 5 in the NIR arm per pair of 
nodded exposures. For SDSS J1212-0123 we achieved a S/N 
of ~ 20 in the UVB arm per exposure, ~ 30 in the VIS arm 
per exposure and ~ 20 in the NIR arm per pair of nodded 
exposures. 



or egress and could not measure the radial velocity ampli- 
tude of the white dwarf. They found that SDSS J1212-0123 
contained a relatively hot (17, 700 ± 300K) low mass white 
dwarf with an active M4 main-sequence companion in an 
8 h 3 m period. 



2.4 Flux calibration 

The ULTRACAM and SOFI light curves were used to flux 
calibrate the X-shooter spectra. We fitted a model to each 
of the light curves (see Section \3.1. 8 1 in order to reproduce 
the light curve as closely as possible. The model was then 
used to predict the flux at the times of each of the X-shooter 
observations (GK Vir shows no stochastic variations or flar- 
ing, whilst small flares were seen in one ULTRACAM obser- 
vation and two SOFI observations of SDSS J1212-0123 but 
were removed before fitting). 

We then derived synthetic fluxes from the spectra for 
the ULTRACAM u , g', r' , i' and z filters as well as the 
SOFI J, H and K filters. We extrapolated the light curve 
models to those bands not covered by our photometry. We 
then calculated the difference between the model and syn- 
thetic fluxes and fitted a second-order polynomial to them. 
This correction was then applied to each spectrum. This cor- 
rects for variable extinction across the wavelength range, as 
well as variations in seeing. 



3 RESULTS 

3.1 SDSS J1212-0123 

SDSS J1212-0123 was initially listed as a q uasar candidate 
from t h e Sloan Digital Sky S urvey (SDSS) bv lRichards et al.l 
(|2004 ). ISilvestri et al.1 ((2006) reclassified it as a white dwarf 
plus main-sequence binary a nd ec lipses were discovered by 
iNebot Gomez-Moran et al.l ([2009) who derived the basic 
system parameters. However, their analysis was limited by 
the fact that they did not resolve the white dwarf ingress 



3.1.1 Spectral features 

Figure [T] shows an average spectrum of SDSS J1212-0123. 
The white dwarf dominates the spectrum at wavelengths 
shorter than 0.55 microns whilst at longer wavelengths the 
spectral features of the secondary star dominate. There are 
also numerous emission lines throughout the spectrum orig- 
inating from the secondary star. Additionally, both Call 
3934A and Mgn 448lA absorption from the white dwarf are 
seen, likely the result of low level accretion from the wind 
of the secondary star. The fact that absorption features are 
seen from both stars allows us to measure the radial veloci- 
ties for both components of the system directly. 



3.1.2 Atmospheric parameters of the white dwarf 

INebot Gomez-Moran et al.l |2009t) determined the temper- 
ature of the white dwarf in SDSS J1212-0123 by decom- 
posing the SDSS spectrum. We computed the average X- 
shooter spectrum of the white dwarf by removing the M- 
dwarf contribution from the individual spectra using obser- 
vations of the M4 star GJ 447, and shifting the residual spec- 
tra to the rest frame of the white dwarf, however this does 
not remove the emission components. We then fit the av- 
er age white dwarf spectrum f ollowing the method outlined 
in iRebassa-Mansergas et al.l (|2007l . l2010f ) . using a grid of 
pure hydroge n model atmos pheres calculated with the code 
described by iKoested l|2010h . We down-weight the cores of 
the Balmer lines since they are contaminated by emission 
from the secondary star, we do not use the Ha? line since 
the secondary star dominates at this wavelength, all other 
lines were used. We find a temperature of 17, 707 ± 35K, 
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Hi 


n 







Table 3. White dwarf absorption features in SDSS J1212-0123. 



3920 3930 3940 4470 4480 4490 

Wavelength [AJ 

Figure 2. The X-shooter spectra of SDSS J1212-0123, averaged 
in the white dwarf rest-frame and normalised (black) along with 
the best-fit white dwarf model (T c ff = 17707, logg = 7.51) and 
log[Mg/H] = —5.8. The CaK line is contaminated by the emission 
from the secondary star (which appears shifted/smeared in the 
rest frame of the white dwarf), however, adopting the same abun- 
dance for Ca provides a reasonable match to the photospheric 
absorption line. 



a surface gravity of log g — 7.51 ± 0.01 and a distance of 
228 ± 5pc, all of which is consist ent with the results from 
iNebot Gomez-Moran et alj (|2009h (and from our light curve 
fit, see Section 2|. Note that these are purely statistical un- 
certainties. 

Next, we fitted the equivalent width of the observed 
Mg II 4481 A absorption line (90 ± 20 mA) by varying the Mg 
abundance in the synthetic spectrum, keeping T e s and log g 
fixed to the values determined above. The best-fit abun- 
dance by number is log[Mg/H] = —5.8 ± 0.1, correspond- 
ing to ~ 4% of the solar value. Figure [2] shows the fit to 
the Mgll448lA line. The CaK line is significantly contam- 
inated by the emission line of the companion star, however, 
Figure [2] illustrates that adopting the Mg abundance also 
for Ca matches the observed CaK line reasonably well. 

Since there is no convection zone acting as a reservoir 
for the accreted elements, the settling times vary throughout 
the atmosphere. At optical depth ~ 2/3, representative for 
the visible spectrum, the timescale for Mg is a few tens of 
days. The diffusion time scale for the te mperature and sur- 
face g ravity derived above is ~ 2 months l|Koester fc Wilkenl 
2006), it is hence plausible to assume accretion-diffusion 
equilibrium. Hence, Xm s Pvm s = constant, with Xm s and 
WM g the mass fraction and the diffusion velocity (relative 
to hydrogen) of Mg, and p the mass density. We com- 
pute at tr oss w 2/3 p = 5.23 x 10 _8 gcm -3 , Xjytg = 
3.76 xlO -5 , and VMg = 8.48 x 10 cms 1 , i.e. a mass flux of 
1.67 x 10 _13 gcm _2 s _1 . Multiplying by the white dwarf sur- 
face area, with Rwd from Table [9j the total mass flux of Mg 
is 2.9 x 10 6 gs _1 - which is equal to the mass accretion rate 
of Mg at the surface of the white dwarf. Assuming that the 
donor star transfers material of solar composition, the total 
accretion rate is then obtained by dividing the Mg rate by 
the mass fraction of Mg in the Sun, i.e. M ~ 4.0 x 10 9 gs _1 
or 6.4 x 10 _17 MQyr _1 . This value is right in the middle 
of the accretion rates calculated for the pre-CVs RRCae, 
UZ Sex, EGUM a LTT560 and S PSS J 12101 0.1+334722.9 
by iDebesI (|2006T l. iTappert et ail l|201ll ), and iPvrzas et all 
l|201ll ), 9 x 1O" 19 M yr _1 to 5 x 10" 15 M r 



Line 



7WD 





(km 


s 




(km s 


- 1 ) 


Can 3933.663 


102.71 


± 


1.32 


37.36 ± 


1.12 


H<5 4101.735 


106.22 


± 


2.21 


35.14 ± 


1.85 


H7 4340.465 


104.28 


± 


1.24 


37.43 ± 


1.08 


Mgn 4481.126 


104.53 


± 


0.78 


35.75 ± 


0.67 


H/3 4861.327 


105.49 


± 


0.93 


35.25 ± 


0.42 


Ha 6562.760 


103.32 


± 


1.13 


38.66 ± 


0.97 



3.1.3 White dwarf radial velocity 

The orbital phases of the X-shooter spectra of SDSS J1212- 
0123 were determined using the ephemeris derived in Sec- 
tion 13.1.61 As previously mentioned both Ca 11 and Mg 11 
absorption from the white dwarf are present. Additionally, 
the cores of the longer wavelength hydrogen absorption lines 
(US to Ha) are narrow and suitable for radial velocity mea- 
surements. 

We measured the radial velocities of the absorption lines 
by simultaneously fitting all of the spectra. We used a com- 
bination of a straight line and Gaussians for each spectrum 
(including a broad Gaussian component to account for the 
wings of the absorption in the case of the Balmer lines) and 
allowed the position of the Gaussians to change velocity ac- 
cording to 

V = 7 + Ksin(2n<l>), 

for each star, where 7 is the velocity offset of the line from its 
rest wavelength and <f> is the orbital phase of the spectrum. 

The parameters determined from the fits to the white 
dwarf absorption features are listed in Table [3] The fits to 
several of the absorption features are also shown in Fig- 
ure [3] the fit to the Ha line is also shown in Figure [4] 
Taking a weighted average of the radial velocities gives a 
radial velocity amplitude of the white dwarf of Kwo = 

104.4 ± 0.5 kms" 1 . 



3.1.4 Secondary star radial velocity 

There are both absorption and emission features originating 
from the secondary star seen in the X-shooter spectra. How- 
ever, the emission features are due to irradiation from the 
white dwarf hence they do not track the true radial velocity 
of the secondary star but the centre of light of the emission 
region, which will be offset towards the white dwarf. There- 
fore we can only directly measure the velocity amplitude for 
the secondary star via the absorption lines. 

Sodium and Potassium absorption lines are present in 
the spectra which we fitted in the same way as the white 
dwarf absorption features. The result of these fits are given 
in Table 3] There were also a number of molecular absorp- 
tion features which we do not fit due to the uncertainty in 
their rest wavelengths and their broad, non-Gaussian pro- 
files. Other atomic absorption features are seen at longer 
wavelengths (for example: the Nal doublet at 2.2 microns) 
but the S/N of these features are too low to reliably fit 
them. A weighted average of the radial velocities of the 
lines in Table [4] gives a radial velocity amplitude of the 
secondary star of K BCC = 168.3 ± 0.3 km s -1 . This is some- 
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Figure 3. Trailed spectra of several lines in SDSS J1212-0123. The grey-scale runs from white (75 per cent of the continuum level) to 
black (125 per cent of the continuum level). The Call 3934A line shows an absorption component from the white dwarf and an emission 
component from the M star. The MgH 448lA line is from the white dwarf, a weak Fel emission line is also visible. Ha shows both 
absorption from the white dwarf and emission from the M dwarf. The K I absorption line and the Na I absorption doublet originate from 
the M star, a telluric correction was applied but artifacts still remain in the Nal trail. The red lines (online version only) show the best 
fits to the lines. 
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Table 4. Secondary star atomic absorption features in 
SDSS J1212-0123. There are several additional atomic absorption 
features in the NIR arm spectra but the S/N of these lines are 
too low to reliably fit them. 



Line 



Ksec 

(kms" 1 ) 



O'sec 

(kms -1 ) 



0.4 0.6 

Orbital Phase 



Ki 7664.899 
Ki 7698.964 
Nal 8183.256 
Nai 8194.824 
Nai 11381.45 
Nai 11403.78 
Ki 11690.219 
Ki 12432.274 
Ki 12522.141 



169.36 
169.04 
168.77 
167.42 
165.24 
166.87 
170.22 
168.63 
167.40 



± 1.13 
±0.83 
±0.45 
±0.48 
±3.86 
± 1.20 
± 1.37 
± 1.01 
± 1.03 



18.27 
19.28 
18.44 
19.52 
21.51 
19.14 
20.86 
19.47 
21.48 



±0.97 
±0.71 
±0.38 
±0.41 
±3.21 
± 1.07 
± 1.24 
±0.88 
±0.72 



Figure 4. Radial velocity fits to the Ho absorption from the 
white dwarf and Nal 8183A absorption from the secondary star 
in SDSS J1212-0123 with residuals plotted below. 



wh at lower than the value of fi " sec = 181 ± 3 km s 1 found 
by iNebot Gomez-Moran et"aT. | (|2009l ). To try and resolve 
this difference we re-fitted the r adial velocity data from 
INebot Gomez-Moran et alJ \200$ ) and found a \ = 25 for 
10 points with 2 variables, implying that they slightly un- 
derestimated their error. Accounting for this, their value is 
consistent with ours to within 2.5<r but, based on our higher 
resolution and additional clean features (not affected by tel- 
luric absorption), we favour our value for -Ksec- We find no 



evidence that the absorption lines are affected by irradia- 
tion from the white dwarf. The equivalent widths of these 
lines do not vary with phase and no effects are visible in the 
radial velocity curves. 

We also measured the radial velocity amplitudes of all 
of the identified emission lines. However, for these lines we 
allow the height of the emission line to vary with orbital 
phase according to 

H = Ho- Hi cos(2tt0), 

which allows the height to peak at phase 0.5, where the 
irradiation effect is largest. This approach gives better fits 
than keeping the height at a fixed value. 

The results of these fits are given in Table \M] in the 
appendix. We find that the emission lines give a range of 
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Table 5. Eclipse times for SDSS J1212-0123. (1) 
iNebot Gomez-Moran et al.l J2009l) , (2) This paper. 



Cycle 
No. 



MJD(BTDB) 
(mid-eclipse) 



Reference 






54104.2092(21) 


(1) 


122 


54145.1854(8) 


(1) 


125 


54146.1929(8) 


(1) 


205 


54173.0628(10) 


(1) 


410 


54241.9170(21) 


(1) 


1455 


54592.9008(14) 


(1) 


3593 


55310.9934268(59) 


(2) 



radial vel ocities, likely due to t he different optical depths of 
the lines (|Parsons et alj l2010al ) . Figure [3] shows the fits to 
the Nal doublet as well as several emission lines. The fit to 
the Nal 8183A line is shown in Figure 3] 



3.1.5 White dwarf's gravitational redshift 

General relativity tells us that the gravitational redshift of 
a white dwarf is given by 



V x = O.635(M/M )(J? / J R)kms- 1 



(1) 



where M and R are the mass and radius of the white dwarf. 
Furthermore, if we know the radial velocity amplitudes of 
the two stars then Kepler's third law tells us 



Mwd = 



PK sec (K W v + K s , 



(2) 



2nG sin 3 i 

where P is the orbital period, and i is the orbital inclination. 
Therefore for a given inclination we can calculate the mass of 
the white dwarf via Eq (J5J) and the radius of the white dwarf 
via a model fitted to the primary eclipse, and thus predict a 
redshift. Hence we can use the measurement of the gravita- 
tional redshift to constrain the inclination by rejecting light 
curve models which do not satisfy this constraint. 

The gravitational redshift can be measured from the 
difference in the velocities of the two components (7wd — 
7scc). 7scc is tightly constrained due to the large number of 
emission lines as well as many absorption features. Taking 
an inverse variance weighted mean of the secondary star line 
velocities (from Tables [4l and lA"Tj) we found 7 SCC = 19.93 ± 
0.06 kms -1 . We calculated 7wd in the same way using the 
values in Table which gave 7wd = 36.0 ± 0.3 km s -1 . We 
found no evidence of pressure shifts in these lines, although 
these are expected t o be small for calcium and magnesium 
l|Vennes et alj|2011af ). 

Using these measurements we determined the gravi- 
tational redshift of the white dwarf to be V z = 16.1 ± 
0.3kms _1 . The true gravitational redshift ofthe white dwarf 
will actually be slightly higher than this value since the mea- 
sured value includes the effects of the secondary star (see 
Section [3. 1.81 for details of these corrections). 



3.1.6 Eclipse time 

We recorded one eclipse of SDSS J1212-0123 with ULTRA- 
CAM. This provides the first high-precision eclipse time for 
this system. The SOFI eclipses are not suitable for long term 
period studies as they are not precisely timed. The eclipse 



time, as well as all previous eclipse times, are listed in Ta- 
ble \5\ The new time is consisten t with the ephemeris of 
INebot Gomez-Moran et al.l |2009) but since it is of higher 
precision we use it to update the ephemeris to 



MJD(BTDB) 



54104.209 17(48) + 0.335 870 93(13)£, 



where we have used barycentric dynamical time (TDB) cor- 
rected for the light travel time to the barycentre of the solar 
system (BTDB). Since this is the first precise eclipse time, 
no long-term period trend is yet visible in the data. 

3.1.7 Spectral type ofthe secondary star 

From their spect ral decomposition 

INebot Gomez-Moran et' al. (2009) determined the spectral 
type of the secondary star in SDSS J1212-0123 to be 
M4 ± 1. We detect the secondary star during the eclipse in 
all bands except the u band giving us multi-colour infor- 
mation. We measure magnitudes for the secondary star of 
g' = 19.73±0.04, i' = 17.384±0.004 and J = 14.949±0.001. 
The i! — J = 2.435 ± 0.004 c olour implies a spe ctral type 
of M4 for the secondary star (Hawlcy ct al. 2002' , Table 3), 
consis tent with the result from Nebot Gomez-Moran et all 
(2009). The distance determined from the spectroscopic fit 
(228 ± 5 pc) means that the absolute J band magnitude of 
the secondary star is Mj = 8 . 2. Usi ng the Mj-spectral type 
relation from lHawlev et all (|2002h . also gives a spectral 
type of M4. 



3.1.8 Modelling the light curves 

The light curves of SDSS J1212-0123 are characterised by 
a deep eclipse of the white dwarf (which gets shallower in 
the longer wavelength bands). There is little variation out 
of eclipse and we found no evidence of a secondary eclipse 
in any of the light curves. 

We fitted all of our photometry (for both SDSS J1212- 
0123 and GKVir) using a code written to produce models 
for the general case of binar ies containing a white dwarf 
(see ICopperwheat et a.1.1 120101 for details). It has been used 
in the study of other white dwarf-main sequence binaries 
l|Pvrzas et al.ll2009l ; |Parsons et al.ll2010al ). The program sub- 
divides each star into small elements with a geometry fixed 
by its radius as measured along the direction of centres to- 
wards the other star. Roche geometry distortion and irra- 
diation of the secondary star are included, the irradiation 
is approximated by crT'f ec = o"T a t c + AF- 1VV where T' scc is 
the modified temperature and T scc is the temperature of the 
unirradiated main-sequence star, a is the Stefan-Boltzmann 
constant, A is the fraction of the irradiating flux from the 
white dwarf absorbed by the secondary star and Fi It is the 
irradiating flux, accounting for the angle of incidence and 
distance from the white dwarf. 

The parameters needed to define the model were: the 
mass ratio, q — M 8CC /Mwd, the inclination, i, the sum of the 
unprojected stellar orbital speeds V a = (Kwd+K ecc ) / sin(i), 
the stellar radii scaled by the orbital separation R scc /a and 
Rwo/a, the unirradiated temperatures, Teff,WD and T e ff, aec , 
quadratic limb darkening coefficients for the both stars, the 
time of mid eclipse, To and the period, P. Note that the 
temperatures are really just flux scaling parameters and only 
approximately correspond to the actual temperatures. 
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Figure 5. ULTRACAM u' , g' and i' and SOFI J band primary eclipses of SDSS J1212-0123 with model fits and residuals. 



Table 6. Parameters from Markov chain Monte Carlo minimisation for SDSS J1212-0123, some fitted, some fixed a priori (those without 
quoted uncertainties), a and b are the quadratic limb darkening coefficients. A is the fraction of the irradiating flux from the white dwarf 
absorbed by the secondary star. 



Parameter 


u' 


a' 


i' 


J 


i (deg) 


86.1 ± 2.2 


85.7 ±0.6 


85.8 ±0.8 


85.5 ±0.8 


r WD /a 


0.0092 ± 0.0002 


0.0092 ± 0.0002 


0.0093 ± 0.0002 


0.0092 ± 0.0002 


^*sec / Q> 


0.171 ±0.016 


0.171 ±0.005 


0.170 ±0.006 


0.172 ± 0.006 


Teff.scc (K) 


2618 ± 112 


2947 ± 26 


3009 ± 41 


3342 ± 42 


"WD 


0.2444 


0.1340 


0.1071 


0.0639 




0.2256 


0.2899 


0.1891 


0.1342 


ttsec 


0.5866 


0.6720 


0.4193 


0.0254 


^scc 


0.2959 


0.2660 


0.4109 


0.4826 


A 


1.60 ±0.50 


0.81 ±0.09 


0.45 ± 0.03 


0.39 ± 0.02 



The light curves are only weakly dependent upon q and 
V s ■ However, we can use them, as well as our measurements 
of the radial velocity amplitudes and the gravitational red- 
shift of the white dwarf, to help constrain the orbital incli- 
nation. This is done by computing the radial velocity am- 
plitudes via 

ifwD = {q/{l + q))V s smi (3) 

K scc = (1/(1 + q))V B sin i, (4) 

we can also compute the masses using 

M = PV^/2-kG, (5) 

where M is the total system mass. The individual masses 
are then 

M WD = M/(l + q) (6) 
M scc = gMwD- (7) 
The orbital separation a is then calculated from 
a = V 3 /2nP, (8) 
allowing us to calculate the radii of the two stars. Combining 



all these calculations yields the masses and radii of both 
stars. We can then use these to calculate the gravitational 
redshifts using 

V z wd = 0.635 (^£ + ^ + (Wsin») 2 

V itwD a J 2c 
T4, scc = o. 635 (^ + ^) + (Wsin^ 

then the value V z = 14, wd — V Z)Sec is equivalent to the mea- 
sured redshift from the X-shooter spectra, which takes into 
account the effects of the secondary star on the measurement 
of the white dwarf's gravitational redshift. Therefore after 
we have generated a model we can compute JfwD, K sec and 
V z from the fitted parameters and either reject or accept 
that model based on how close they are to the measured 
values. 

We used the Markov Chain Monte Carlo (MCMC) 
method t o determine the distributions of our model pa- 
rameters (|Press et all 120071 ). The MCMC method involves 
making random jumps in the model parameters, with new 
models being accepted or rejected according to their proba- 
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Wavelength (microns) 

Figure 6. Averaged X-shooter spectrum of GKVir. The SDSS u'g'r'i'z' and UKIDSS JH magnitudes and filter widths are also shown, 
the UKIDSS J band observations were made during the primary eclipse. A zoom in on the white dwarf features are shown inset with 
the narrow Hen absorption feature labelled. 



bility computed as a Bayesian posterior probability. In this 
instance this probability is driven by a combination of x 2 
and the prior probability f rom our spectroscopic co nstraints 
(see IParsons et al.l l2011bl and IPvrzas et al.l l201ll for more 
details of the MCMC fitting process) . 

For fitting the light curves of SDSS J1212-0123, we 
phase folded the data and kept the period fixed as 
one. We also kept the temperature of the white dwarf 
fixed at 17,900K (see Section I3.1.2|l . For the secondary 
star we used quadratic limb darkening coefficients from 
IClaret fc Bloemenl <|201ll ) for a T cff = 3000, log g = 5 main 
sequence star. For the white dwarf we calculated quadratic 
limb darkening coefficients from a white dwarf model with 
Twd = 17,900 and log g — 7.53 based on our spectroscopic 
fits, folded through the ULTRACAM u, g' , i' and SOFI J 
filter profiles. For both stars we quote the coefficients a and 
b where I(n)/I(l) = 1 — o(l — (J-) — b(l — n) 2 , where fi is the 
cosine of the angle between the line of sight and the surface 
normal. We kept all limb darkening parameters fixed. 

Table[H]lists the best fit parameters from the light curves 
and their la uncertainties, we also list the limb darkening 
coefficients used for each band. The results from all four 
bands are consistent. Figure shows the fits to the primary 
eclipses in each band and the residuals to the fits. 



3.2 GKVir 



GKV ir (PG 1413 + 015) was discovered by iGreen et al.l 
(|l978h fr om the Palomar-G reen survey for ultraviolet-excess 
objects l|Green et al.lll986T ). An eclipse wa s recorded dur- 
ing a subsequent spectroscopic observat ion. Fulbright et all 
(1993) combined the photometry from IGreen et all (I1978T I 
and high resolution spectroscopy to constrain the system 
parameters. They found that GKVir contains a hot DAO 
white dwarf with a low-mass M3-5 main sequence compan- 
ion in an 8 h 16 m period. However, their lack of radial velocity 
information limited their analysis. 



3. 2. 1 Spectral features 

Figure HJ] shows an average X-shooter spectrum of GK Vir. 
The hot white dwarf dominates the spectrum at wavelengths 
shorter than 1 micron. In the J band the secondary star 
contributes roughly 50 per cent of the overall flux. Hy- 
drogen Balmer absorption lines from the white dwarf are 
visible as well as narrow Hell 4686A absorption from the 
wh ite dwarf, making it a DAO white dwarf, as mentioned 
by iFulbright et al.l l|l993l ) . Numerous emission lines origi- 
nating from the heated face of the secondary star are seen 
throughout the spectrum. Nal (8183A, 8195A) absorption 
originating from the secondary star is seen before phase 0.25 
but decreases in strength towards phase 0.5 due to the in- 
creased ionisation and therefore cannot be used to measure 
the radial velocity amplitude of the secondary star. No other 
absorption features for the secondary star are visible. 



3.2.2 Atmospheric parameters of the white dwarf 

We fitted the X-Shooter spectrum, averaged in the white 
dwarf rest-frame, to obtain an estimate of the effective tem- 
perature. As for SDSSJ1212-0123 (see Section l3~i~2)l . we 
fitted the normalised Balmer lines, including H/3 to H£ and 
down-weighting the regions that are noticeably contami- 
nated by emission lines. The best-fit gives a temperature 
of 55995 ± 673K, a surface gravity of log g — 7.68 ± 0.04 and 
a distance of 550 ± 20 pc. We note that the quoted errors 
are purely statistical, and that systematic effects are very 
likely affecting these results. The X-Shooter data were ob- 
tained near the quadrature phases (r/> ~ 0.25 and 0.75, in 
order to measure Kwu and K scc ) , which results in a signifi- 
cant contamination by emission from the strongly irradiated 
inner hemisphere of the secondary star. Varying the wave- 
length range around the lines that is down-weighted affects 
the resulting temperature by several 1000 K. For compari- 
son, we also fitted the lower-resolution SDSS spectrum of 
GKVir, obtained near superior conjunction of the compan- 



© 2011 RAS, MNRAS 000,ITJfT7l 



10 S. G. Parsons et al. 



Ca I 
3934 




He 

587£ 



Ha 
6563 



-400 -200 200 400 - +0 ° - 200 



400 -400 -200 



400 -400 -200 



400 -400 -200 



Velocity (km/s) 

Figure 7. Trailed spectra of several lines in GK Vir. The grey-scale runs from white (75 per cent of the continuum level) to black (125 
per cent of the continuum level). The Can 3934 A line shows emission from the M star as well as weak interstellar absorption which 
shows no radial velocity variations. The Hell 4686A line originates from the white dwarf. The H/ j and Ha emission from the M star 
shows inverted cores and asymmetric profiles similar to those seen in NNSer jParsons et aljfeoiOal) . The red lines (online version only) 
show the best fits to the lines. 



ion, and find T off = 52258 ± 3131 K, log g = 7.66 ± 0.18, and 
d = 509 ± 61 pc. 

iFulbright etaH l| 19931 ) analysed a blue spectrum of 
GKVir obtained near the eclipse ((/> ~ 0.02), i.e. when the 
heated inner hemisphere of the companion contributes least 
to the observed flux. By fitting the H/3-H£ Balmer lines and 
the Hen 4686A line they determined T cff = 48800 ± 1200 K 
and logp = 7.70 ± 0.11. They also modelled the single 
available far-ultraviolet spectrum of GKVir, and obtained 
TL50000K. A major limitation of this spectrum obtained 
with IUE was that the photospheric Lya line was nearly 
completely filled in by geocoronal emission. 

We conclude that the effective temperature and dis- 
tance of the white dwarf GKVir remains somewhat uncer- 
tain, T cff ~ 50000 K and d ~ 500 ± 50 pc. However, it is 
reassuring that all the spectroscopic measurements of the 
surface gravity are consistent with the value that is deter- 
mined from the light curve fit (see Section [4j> . 

Fitting the Hen 4686A absorption line from the white 
dwarf in GKVir gives a helium abundance of log[He/H] = 
—2.8 ± 0.3 by numbers. Assuming that the secondary star 
transfers material of solar composition, the total accretion 
rate is M ~ 1.4 x 10 9 gs _1 or 2.2 x 10" 17 M Q yr" 1 , once 
again consistent with other pre-CVs. 
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Figure 8. Radial velocity fits to the Hell 4686A absorption from 
the white dwarf and H/3 emission from the secondary star in 
GKVir with residuals plotted below. The emission component 
does not track the centre of mass of the secondary star (see Sec- 
tion ESS}- 



3. 2. 3 White dwarf radial velocity 

Due to the high temperature of the white dwarf, the Balmer 
lines lack narrow cores and are therefore unsuitable for ra- 
dial velocity work. We measured the radial velocity of the 
white dwarf in GKVir from the Hen 4686A absorption 
line using the same technique as we used for SDSSJ1212- 
0123. Figure [7| shows a trail of the Hen line and Fig- 



ure [8] shows the fitted radial velocity curve, which gives a 
value of -Kwd = 38.6 ± 0.8 km s" 1 and a velocity offset of 
7wd = -27.2 ± 0.7 km s" 1 . 



3.2.4 Emission lines 

A large number of emission lines are seen in the spectrum of 
GKVir caused by the heating of the secondary star by the 
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Figure 9. Left: radial velocities of the emission lines from the secondary star in SDSS J1212-0123. The measured value of centre of mass 
radial velocity (K Bec ) is also shown. The two dashed lines denote the radial velocity that an optically thick and optically thin line would 
have. The solid line at 133 kms -1 is a weighted average of the radial velocity amplitudes of the Fel lines, these appear to be optically 
thin. The radial velocity amplitudes of the hydrogen Balmer lines and the calcium H and K lines are far above the optically thick limit 
implying that these lines originate, at least in part, from stellar activity. Right: a similar plot but for GKVir. However, in this case we 
were unable to directly measure K Bec from the spectroscopy, but by assuming that the Fel lines are optically thin we get a consistent 
solution. The much tighter limits and smaller correction in GK Vir imply that the secondary star is quite small. 



white dwarf. However, the secondary star in GKVir receives 
a much larger irradiating flux to that in SDSS J1212-0123 
due to the higher temperature of the white dwarf, therefore 
a larger number of lines are present and from higher ionised 
states. Unfortunately the strength of the Nal absorption 
doublet is too low after phase 0.25 and, combined with the 
lack of any other absorption features from the secondary 
star, means that we are unable to get a direct measurement 
of the radial velocity amplitude of the centre of mass of the 
secondary star. We determine the radial velocities and off- 
sets of all of the emission lines identified using the same 
method used for SDSS J1212-0123, these are listed in Ta- 
ble lA2l in the Appendix. We use these results to determine a 
velocity offset of 7 ae c = —47.35 ± 0.05 km s -1 , giving a mea- 
sured redshift for the white dwarf of V z = 20.2 ± 0.7 kms -1 . 



3.2.5 K scc correction 

As previously noted the emission lines in the X-shooter spec- 
tra of GK Vir cannot be used to directly measure K sec , 
needed for accurate mass determinations. We need to de- 
termine the deviation between the reprocessed light centre 
and the centre of mass for the secondary star. The radial 
velocity of the centre of mass (K scc ) is related to that of the 
emission lines (Kemis) by 

Ksec = i-/(i7kc/« 

ijParsons et al.ll2011bl ). where / is a constant between and 
1 which depends upon the location of the centre of light. 
For / = the emission is spread uniformly across the entire 
surface of the secondary star and therefore the centre of light 
is the same as the centre of mass. For / = 1 all of the flux 
is assumed to come from the point on the secondary star's 
surface closest to the white dwarf (the substellar point). 

The centre of light for an emissi on line is related to 
the optical depth of the emission l|Parsons et al.l l2010al ; 



iParsons et ai]|2011bl ). Optically thick emission tends to be 
preferentially radiated perpendicular to the stellar surface, 
therefore at the quadrature phases, we will see the limb of 
the irradiated region more prominently (compared to the re- 
gion of maximum irradiation) than we would otherwise. This 
will lead to a higher observed semi-amplitude and hence a 
smaller correction factor is needed. The reverse is true for op- 
tically thin lines where the emission is radiated equally in all 
directions, hence emission from the substellar point becomes 
enhanced at quadrature, leading to a low semi-amplitude 
and a larger correction factor. The correction factor for an 
optically thin line (assuming emissivity proportional to the 
incident f lux) is / = 0.77 and / = 0.5 for optically thick 
emission l|Parsons et al.ll2010al ). 

We can estimate the optical depths of the emission lines 
in GK Vir using our observation of SDSS J1212-0123 and as- 
suming a similar behaviour in the lines. We can determine 
the optical depth of the lines in SDSS J1212-0123 because 
in this case we have a direct measurement of K scc . There- 
fore, we can reverse Eq (|lip to determine the optical depths. 
The left hand panel of Figure [9] shows where the emission 
lines lie with respect to the optically thin lower limit and 
the optically thick upper limit. We find that the measured 
radial velocities of the Fel lines are consistent with them 
being optically thin. The majority of the emission lines in 
SDSS J1212-0123 appear to lie somewhere between optically 
thick and optically thin with the exception of the hydrogen 
Balmer lines and the calcium H and K lines. For these lines 
the measured radial velocity is far higher than the optically 
thick limit, implying that the emission is more uniformly 
spread over the surface of the secondary star, this is most 
likely caused by stellar activity. This is consistent with the 
observations of several small flares during our photometry 
and implies that the secondary star in SDSS J1212-0123 is 
an active star. 

From our analysis of the emission lines in SDSS J1212- 
0123 we make the assumption that the Fel lines in GKVir 
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Table 7. Eclipse tim es for GKVir. (1) iGreen et al] [|l978h . (2) 
iParsons et"afl feOlObh . (3) this paper. 



3004 



2006 



2008 



2010 



Cycle 


MJD(BTDB) 


Reference 


No. 


(mid-eclipse) 




-67 


42520.26747(1) 


(1) 


-32 


42532.31905(2) 


(1) 


-29 


42533.35204(9) 


(1) 





42543.33769(1) 


(1) 


3 


42544.37068(1) 


(1) 


851 


42836.36314(6) 


(1) 


1966 


43220.29202(12) 


(1) 


2132 


43277.45101(6) 


(1) 


2896 


43540.51972(12) 


(1) 


28666 


52413.9255716(9) 


(2) 


29735 


52782.0152272(9) 


(2) 


29738 


52783.0482185(7) 


(2) 


30746 


53130.1336878(27) 


(2) 


32706 


53805.0221154(23) 


(2) 


32709 


53806.0551129(12) 


(2) 


34054 


54269.1800868(3) 


(2) 


37069 


55307.3375852(11) 


(3) 



are also optically thin. We can then use their radial velocities 
to predict K Bec using Eq pip . We adopt a value for the 
radial velocity amplitude of an optically thin line of Ktmu = 
204 ± 2kms _1 . We use this result as a prior constraint in 
our light curve fitting. The right hand panel of Figure [§] 
shows the radial velocity amplitudes of the emission lines in 
GK Vir, the dashed lines are based on the results of our light 
curve fitting (see section 2. 8 1) . The spread in the emission 
line radial velocities in GKVir is relatively small compared 
to SDSS J1212-0123 because of the small relative size of the 
secondary star (R BCC /a). 



3.2.6 Eclipse times 

We recorded one new high precision eclipse time for GK Vir 
which is listed in Table along with all previous eclipse 
times. We update the ephemeris to 



MJD(BTDB) 



42543.337 9121(33) 

+ 0.344 330 832 742(99)£, 



which is consistent with previous studies. Figure [10] shows 
the difference between the observed eclipse time and the 
calculated eclipse time based on our new ephemeris. Our 
new eclipse time shows a clear deviation from linearity, al- 
though the magnitude of t he period change is s mall com- 
pared to other systems (see IParsons et al.ll2010bh . The sec- 
ondary star in GK Vir is able to drive this small period 
change (0.00124 seconds in ~ 8 years) via Applegate's mech- 
anism (Apple gatelll992l ). although this is only just the case 
if we used t he modified version of A pplegate's mechanism 
presented by iBrinkworth et al] (|2006t ) which takes account 
of the role of the inner part of the star in counterbalancing 
the outer shell. A third body in orbit around the system 
may also be the cause of the period change, as h as recently 
been proposed for NN Ser (|Beuermann et alboiOl ). however, 
longer term monitoring is required in order to discover the 
true cause of this period change. 



o 
o 



3x10* 



Cycles 



3.5x10* 



Figure 10. Observed-Calculated (O-C) plot for the eclipse times 
of GK Vir. Our new point shows a departure from linearity in the 
eclipse times. The secondary star is just able to supply the energy 
required to drive this period change via Applegate's mechanism 
however, a third body in orbit around the system is also a possible 
explanation. 



3.2.7 Spectral type of the secondary star 

We detect the secondary star in GKVir in the r', i' and J 
band eclipses. We measure mag nitudes of r = 21.72 ± 0.03. 
i' = 19.98 ± 0.01 and J = 17.59 ± 0.05. The r' - i' = 
1.7 4 ± 0.03 colour is consistent with a spectral type of M4.5 
whil st the i! — J = 2.39 ± 0.05 is closer to a spectral type of 
M4 l|Hawlev et al.ll2002l ). Using the distance from our spec- 
troscopic fit (550±20pc) the secondary star has an absolute 
J band magnitude of Mj = 8.9, giving it a spectral type 
of M4.5 l|Hawlev et alj|2002h , therefore we adopt a spectral 
type of M4.5 ± 0.5. 



3.2.8 Modelling the light curves 

The light curves of GK Vir show a deep eclipse of the white 
dwarf even in the J band. GK Vir shows a small reflection ef- 
fect out of eclipse, however no secondary eclipse is detected. 
We fit the light curves of GKVir in the same way as those 
of SDSS J1212-0123 (see Section |3~1~5|| . However, since we 
lack a direct measurement of K soc , we use our constraint on 
if Thin and Eq (|11[) as well as our .Kwd and gravitational 
redshift measurements to constrain the inclination. 

We phase binned the data, using smaller bins on the 
ingress and egress, and kept the period fixed as one when 
fitting the light curves. We also kept the temperature 
of the white dwarf fixed at 55,995K. For the secondary 
star we used quadratic limb darkening coefficients from 
IClaret fc Bloemenl (|201ll ) for a T cff = 3000K, log g = 5 main 
sequence star. For the white dwarf we calculated quadratic 
limb darkening coefficients from a white dwarf model with 
Twd = 55,995 and logg = 7.68 based on our spectroscopic 
fits, folded through the ULTRACAM u , g', r' and i' and 
SOFI J filter profiles. For both stars we quote the coefficients 
a and b where 7(/x)/7(l) = 1 — a(l — fi) — 6(1 — /i) 2 , where 
H is the cosine of the angle between the line of sight and 
the surface normal, we kept all limb darkening parameters 
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Orbitol Phose 

Figure 11. ULTRACAM u', g' , r' and i' and SOFI J band primary eclipses of GKVir with model fits and residuals. The light curves 
have been binned up by a factor of 5 although smaller bins were used on the ingress and egress features, which leads to the small spikes 
seen in the residuals at the ingress and egress. 



Table 8. Parameters from Markov chain Monte Carlo minimisation for GKVir, some fitted, some fixed a priori (those without quoted 
uncertainties), a and b are the quadratic limb darkening coefficients. A is the fraction of the irradiating flux from the white dwarf absorbed 
by the secondary star, this was kept fixed for the r' band since there is little coverage out-of-eclipse in this band. 



Parameter 


u' 


a' 


r' 


i' 


./ 


i (deg) 


89.5 ±0.7 


89.5 ±0.6 


89.3 ±0.9 


89.2 ±0.9 


88.8 ± 1.4 


r WD /a 


0.0093 ± 0.0002 


0.0094 ± 0.0002 


0.0093 ± 0.0003 


0.0094 ± 0.0003 


0.0097 ±0.0006 


rscc/a 


0.085 ±0.002 


0.085 ±0.002 


0.086 ± 0.002 


0.086 ± 0.003 


0.088 ± 0.006 


T e s, S cc (K) 


3113 ±31 


3270 ± 12 


3224 ± 28 


3508 ± 40 


4117 ± 139 


»WD 


0.0769 


0.0594 


0.0505 


0.0446 


0.0257 


bwo 


0.1393 


0.1165 


0.0896 


0.0736 


0.0490 


ttsec 


0.5866 


0.6720 


0.6364 


0.4193 


0.0254 


frsec 


0.2959 


0.2660 


0.2521 


0.4109 


0.4826 


A 


1.055 ±0.039 


0.551 ±0.006 


0.400 


0.263 ± 0.009 


0.245 ±0.032 




221.0 ±2.1 


221.2 ±2.1 


221.1 ±2.1 


221.3 ±2.1 


222.5 ±2.3 



fixed. For our r' band light curve we lack any out of eclipse 
information (barring that immediately before and after the 
eclipse) hence A, the fraction of irradiating flux absorbed 
by the secondary star, is unconstrained. Therefore, we fix 
this value at 0.4; this parameter has no effect on the radii 
or inclination. 

Figure [11] shows the fits to the light curves and the 
residuals and our final fitted parameters are listed in Table [8] 
We also list the limb darkening coefficients used for each 
band. 



4 DISCUSSION 

Our light curve fits combined with Eqs ©-(JSJl yield direct 
measurements for the masses and radii of both components 
in SDSS J1212-0123 and GKVir. For our final values and 
uncertainties we combine the results of all our light curves, 



however, since they are all constrained by the same spec- 
troscopic information the uncertainties in each light curve 
fit are not independent. Therefore we combine the results 
from each light curve in an optimal way by calculating the 
minimum possible error on each parameter purely from the 
spectroscopic constraints and then combining parameters al- 
lowing for the correlated and random noise components. Our 
final parameters for both systems are listed in Table The 
secondary star's shape in both systems is slightly distorted 
due to the presence of the nearby white dwarf, therefore 
Table |§] lists the radius of the secondary star in various di- 
rections. This is a very minor effect in GKVir whilst the 
effect is somewhat larger for SDSS J1212-0123, though in 
both cases the variations are smaller than the uncertainty 
on the radius from the light curve fits. For our final discus- 
sions we adopt the volume-averaged radii. 

Figure [121 shows the mass-radius plot for white dwarfs. 
The position of the white dwarfs in SDSS J1212-0123 and 
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Table 9. System parameters. The surface gravities quoted arc from the spectroscopic fit. They are consistent with the measured masses 
and radii. The accretion rate is that of the material from the wind of the secondary star onto the white dwarf 



Parameter 



SDSS J1212-0123 



GKVir 



Period (days) 

Inclination 

Binary separation 

Mass ratio 

WD mass 

Sec mass 

WD radius 

Sec radius polar 

Sec radius sub-stellar 

Sec radius backside 

Sec radius side 

Sec radius volume- averaged 

WD log g 

WD temperature 

%D 

Vz,WD 

Sec spectral type 
Distance 
Accretion rate 



0.335 871 14(13) 
85.7° ± 0.5° 
1.815 ±0.003 R 
0.620 ±0.001 
0.439 ± 0.002 M 
0.273 ±0.002 M 
0.0168 ± 0.0003 R Q 
0.304 ± 0.007 R 
0.310 ±0.007 R 
0.309 ± 0.007 R 
0.306 ± 0.007 R 
0.306 ± 0.007 R 
7.51 ± 0.01 
17, 707 ± 35K 
104.4±0.5kms~ 1 
168.3 ±0.3kms~ 1 
16.1 ± O.Skms" 1 
M4 

228 ±5pc 



i.4 x 10" 



'M 



©y r " 



0.344 330 832 742(99) 
89.5° ± 0.6° 
1.82 ±0.01 R 
0.174 ±0.004 
0.564 ±0.014 M 
0.116 ±0.003 M 
0.0170 ± 0.0004 R Q 
0.154 ±0.003 R 
0.156 ±0.003 R 
0.156 ±0.003 R 
0.155 ±0.003 R 
0.155 ±0.003 R 
7.68 ± 0.04 
55, 995 ± 673K 
38.6 ±0.8 km s^ 1 
221.6 ±2.0kms" 1 
20.2 ±0.7^8-! 
M4.5±0.5 
550 ±20 pc 
2.2 X lO-^Mgjyr- 1 




Mass (M Q ) 

Figure 1 2 ■ Mass-radius plot fo r white dwarfs. Black points 
are from IProvencal et al. | il99Si ). IProvencal et al.l i2002T l and 
ICasewell et al. I (I2009T ). The square points are common proper- 
motion systems, the triangular points are visual binaries. White 
dwarfs measureme nts from PCEBs ar e shown in red (onlin e 
version only) from JO'Brien et alj|200ir) , iParsons et"al] ll2010ah . 
iPvrzas et alj 1 120111) and this work. The grey lines are CO core 
models with the temperatures labelled in units of 10 3 K and 
with hydrogen layer thickn esses of Mh/Mwd = 10~ 4 from 
iBenvenuto k Althauj dl999l ') ■ The mass and radius of the white 
dwarf in GK Vir agree well with the 50, 000K model. The blue 
lines (online version only) are He core models with a temperature 
of 18, 000K and varying hydrogen layer thickn esses labelled by th e 
exponent of the hydrogen layer fraction, from lPanei et al. N2007t) . 
The mass and radius of the white dwarf in SDSS J1212-0123 is 
consistent with the CO core models for its temperature but we can 
rule out a CO core on evolutionary grounds. Therefore the white 
dwarf in SDSS J1212-0123 is consistent with the He core models 
only if it has a very thin hydrogen envelope (Mh/Mwd ^ 10 — 6 ). 



GKVir are shown as well as other accurate white dwarf 
mass-radius measurements. The measured mass and ra- 
dius of the white dwarf in GKVir are consistent with a 
carbon-oxygen (CO) core white dwarf of the same tem- 
perature, with a thick hydrogen envelope (Mh/Mwd = 
1CP 4 ). The measured mass and radius of the white dwarf 
in SDSS J1212-0123 are also consistent with a CO core 
white dwarf with the same temperature and a thick hy- 
drogen envelope (Mh/Mwd = 10 -4 ). However, although 
it is possible to create CO core white dwarfs with masses 
< 0.5 MfTi via considerable mass l oss al o ng the red giant 
phase ifPrada Moroni fc Stranierol 120091 ; IWillems fc Kolbl 
l2004l;lHan et al.ll2000l ). doing so in a binary system requires 
a large initial mass ratio and results in a widening of the 
orbital separation, hence we would not expect a CO core 
white dwarf with a mass < 0.5 Mq in a close binary sys- 
tem. Therefore the white d warf in SPSS J12 12-0123 must 
have a He core, as noted bv lShen et all (|2009h . 

Several He core mass-radius relations are shown in Fig- 
ure Q2] for a white dwarf with a temperature of 18, 000K 
and varying hydrogen layer thicknesses. The white dwarf in 
SDSS J1212-0123 is consistent with these relations only if 
it has a very thin hydrogen envelope (Mh/Mwd ^ 10 -6 ). 
Figure [13] shows a zoomed in version of the mass radius 
plot for the white dwarf in SDSS J1212-0123 as well as the 
same He core models as Figure 1121 The black line shows 
the range of possible masses and radii that the white dwarf 
in SDSS J1212-0123 could have based on the radial veloci- 
ties and the primary eclipse shape (i.e. no inclination con- 
straints). The numbers plotted along this line are what the 
measured gravitational redshift of the white dwarf would 
need to be in order to give that mass and radius. Previ- 
ous studies have found that the spectroscopic gravitational 
redshift measurements are usually slightly inc onsistent with 
mass - radius measurements v i a other methods fpyrzas et al.l 
l201ll ; IParsons et all l2010al ; iMaxted et all I2007T I meaning 
that our inclination constraints may be slightly incorrect. 
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Mass (M Q ) 

Figure 13. Mass-radius plot for the white dwarf in SDSS J1212- 
0123. The blue lines (online version only) are He core models 
with a temperature of 18, 000K and varying hydrogen layer thick- 
nesse s labelled by the exponent of the hydrogen layer fraction, 
from iPanei et al.l (120071 ) . The black line shows the possible mass 
and radius range without any constraint on the inclination with 
90° at the top. The numbers along this line correspond to the 
gravitational redshift that the white dwarf would have at that in- 
clination (in kms -1 ). Our measured redshift (16.1 ± O.Skms^ 1 ) 
means that the white dwarf is consistent with models which have a 
thin hydrogen layer. However, even without any inclination con- 
straints it is not possible for the white dwarf to have a thick 
(M H /M WD > 2 X 10 -4 ) hydrogen layer. 



Figure 14. Mass-radius plot for low-mas s stars. Black point s 
are high-prec i sion m easu rements taken fro m lKnigge et al.l l l201ll) , 
ICarter et all d201lh and lOfir et al.l ll201ll). T he solid line is the 
4.0-Gyr isochrone from iBaraffe et al .1 dl998h whilst the dashed 
line is a 4.0-Gyr model from lMoralesetaL 1 120101) which includes 
the effects of magnetic activity. Low-mas s stars from PCEBs ar e 
shown in red (onlin e version only) from iParsons et al.l j2010a| N | , 
Pv rzas et ah! (1201 lh and this work. The blue points (online ver- 
sion only) show the radii of these stars after correcting for irra- 
diation (negligible for SDSS J1210+3347 and SDSS J1212-0123). 
The secondary star in GKVir is consistent with model predic- 
tions after correcting for irradiation, however the secondary star 
in SDSS J1212-0123 is still oversized. 



However, Figure [13] shows that it is not possible for the 
white dwarf to have a thick hydrogen envelope (Mh/Mwd > 
2 x 10~ 4 ) even if the inclination is 90°. Current evolution- 
ary models are unable to create He core white dwarfs with 
such thin hydrogen envelopes, the thi nnest envelopes are 
of the order of Mh/M W d ~3x 10~ 4 l|Althaus et al.ll2009l : 
ISarna et al.1l2000l ;l Driebe et al.|[l99^ ) meaning that, assum- 
ing our measured V z value is accurate, either SDSS J 1212- 
0123 has had a very unusual evolutionary history or that 
current evolutionary models of He core white dwarfs are in- 
complete and overestimate their size. Additional He core 
white dwarf mass-radius measurements should show if this 
is the case. 

Figure [14] show the mass-radius plot for low mass stars. 
The masses and radii of the secondary stars in SDSS J1212- 
0123 and GKVir are marked as well as other precise mea- 
surements. The mass and radius of the secondary star in 
SDSS J1212-0123 show that it is over-inflated for its mass 
compared with evolutionary models by ~ 12 per cent. This 
discrepancy can be reduced to ~ 9 per cent if the sec- 
ondary star is active. In this case the radius may be overes- 
timated due to the effects of polar spots, it also increases to 
compensate for loss of radiative efficiency due to starspots 
and due to a strong magnetic field caused b y rapid rota- 
tion l|Morales et al.ll2010l ; Ichabrier et alj|2007h . We have ev- 
idence from the emission lines and light curves that the sec- 
ondary star in SDSS J1212-0123 is indeed active, therefore 



these effects can explain some of the discrepancy however, 
the star remains oversized. The secondary star in GKVir is 
also oversized by ~ 9 per cent, this drops to ~ 6 per cent 
if it is active, though we have no evidence of activity from 
this star. 

Taking the effects of rotational and tidal perturbations 
fully into account only causes an increase in the radii of 0.1 
per cent for GK Vir a nd 0.4 per cent for SDSS J1212-0123 
(Sirotkin & Kimll2009l ) which is not enough to explain the 
discrepancy. 

Irradiation by the white dwarf can cause the sec- 
ondary star to become inflated by effectively blocking the 
energ y outflow through the surface layers l|Ritter et al.l 
2000). For NNSer correcting this effect brought the mass 
and radius measurements into a greement with evolution- 
ary models (|Parsons et al.l [2010al ). For the secondary star 
in G K Vir we find that irradiation increases its radius by 
5.6% (|Ritter et al.l l2000l : lHameurv fc Ritterlll997l) . enough 
to bring it into agreement with evolutionary models. How- 
ever, for the far less irradiated secondary star in SDSS J1212- 
0123 we find an increase of only 0.4% meaning that it is still 
overinflated. The corrected radii are shown in Figure [14] 



5 CONCLUSIONS 

We have used a combination of ULTRACAM and SOFI 
photometry and X-shooter spectroscopy to measure pre- 
cise masses and radii for both components of the eclips- 
ing PCEBs SDSS J1212-0123 and GKVir. In both cases we 
use measurements of the gravitational redshift of the white 
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dwarf to constrain the orbital inclination. We were able to 
directly measure the radial velocity amplitudes of both stars 
in SDSS J1212-0123 but only the white dwarf in GKVir. 
However, we use our results from SDSS J1212-0123 and in- 
formation from the emission lines to determine the correc- 
tion factor to apply to the emission line radial velocities in 
GK Vir to measure the centre of mass radial velocity of the 
secondary star. The mass and radius of the 0.564 Mq white 
dwarf in GKVir are consistent with evolutionary models 
for a CO core white dwarf with a thick hydrogen envelope. 
The mass and radius of the white dwarf in SDSS J1212- 
0123 place it in the crossover region between a CO core 
white dwarf and a He core white dwarf but we can exclude 
a CO core on evolutionary grounds. This means that the 
white dwarf in SDSS J1212-0123 is the first He core white 
dwarf with precise mass-radius measurements however, it is 
under-inflated for its temperature unless it has a very thin 
(Mh/Mwd ^ 10 -6 ) hydrogen envelope, which evolution- 
ary models are unable to produce. The mass and radius of 
the secondary star in GKVir are consistent with evolution- 
ary models after correcting for the effects of irradiation by 
the white dwarf, however, the secondary star in SDSS J1212- 
0123 is overinflated compared to theoretical predictions even 
after the effects of activity, rapid rotation and irradiation are 
taken into account. 
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Table Al. Secondary star emission lines in SDSS J1212-0123. After measuring the radial velocity of a line, the spectra were shifted to 
remove the motion and averaged. The Equivalent widths were measured from this averaged spectrum. 



Line 




■7sec 

(km/s) 




-Kmc as 

(km/s) 




Equivalent 
Width (mA) 


Mgl 3829.355 


18 


.31 


± 


2 


.52 


141. 


.19 


±3. 


.68 


72 ± 3 


tit T oooo onn 

Mgl 3832.299 


22 


.90 


± 


3 


.20 


135. 


.22 


± 4. 


.12 


60 ± 3 


tv t T oooo onn 

Mgl 3838.292 


17 


.09 


± 


2 


.25 


138. 


.27 


±3. 


.17 


67 ± 3 


T7 1 .. T O O T O ft T O 

lei 3878.573 


18 


.74 


± 


1 


.89 


132. 


.01 


± 2. 


.41 


83 ± 2 


tto ooon n tr c 

Hs 3889.055 


17 


.47 


± 


1 


.73 


157. 


.82 


± 2. 


.18 


74 ± 3 


o; T inner trno 

bi I 3905.523 


22 


.79 


± 





.33 


140. 


.12 


± 0. 


.43 


248 i 2 


tti_ T o n no ni"i 

le I v>yzz.911 


20 


.92 


± 


1 


.86 


132. 


.93 


±2. 


.32 


41 ± 2 


t - . onno / 1 l ; • > 

lei 3928.083 


15 


.52 


± 


2 


.49 


131. 


.28 


± 3. 


.52 


37 ± 2 


/ 1 _ TT onoo cr"~> 

Gall 3933. 6b3 


22 


.31 


± 





.38 


153 


.32 


±0. 


.18 


386 ± 3 


A l T on a a r\c\c 

Al I oy44.00b 


25 


.21 


± 


1 


.41 


141. 


.63 


± 1. 


.99 


48 ± 2 


A l T onci rr on 

Al I o9ol.5z0 


19 


.95 


± 


1 


.01 


140. 


.74 


± 1. 


.35 


72 ± 2 


Ua II oybo.4b9 


21 


.05 


± 





.53 


153. 


.32 


± 0. 


.18 


349 ± 3 


tt onTn ot /i 

He 3970.074 


17 


.03 


± 


2 


.88 


160. 


.10 


±2. 


.13 


84 ± 2 


TT r ^in~i TOP 

rlo 4101. 7oo 


21 


.62 


± 


1 


.23 


157. 


.86 


± 1. 


.53 


262 i 4 


TT1„ T A 1 a n O O 

le I 414o.obo 


22 


.17 


± 


2 


.91 


143. 


.06 


±4. 


.15 


31 ± 2 


T71 T (ill C 10/1 

lei 4zlb.lo4 


20 


.09 


± 


2 


. 30 


137. 


.76 


± 3. 


.16 


27 ± 2 


T /inor too 

Ua I 4zzb.7z8 


22 


.94 


± 


1 


.59 


138. 


.71 


±2. 


.13 


91 ± 2 


tt. , /i o /i n a a r 
H*y 4o40.4bo 


17 


.31 


± 





.60 


160. 


.76 


±0. 


.76 


320 i 3 


TT 1 — T- /i o T tr no/ 1 

lei 4375.980 


21 


.88 


± 


1 


.73 


136. 


.22 


±2. 


.20 


54 ± 2 


T?„ T i 1 o t o ~i n 

le I 44z7.ol0 


19 


.64 


± 


1 


.57 


134. 


.77 


± 2. 


.09 


49 ± 2 


le I 44b 1. boo 


20 


.38 


± 


2 


.23 


141. 


.02 


±2. 


.91 


47 ± 2 


TT 1 .. T A ion 1 TO 

tel 4482.170 


20 


.86 


± 


1 


.09 


129. 


.57 


± 1. 


.22 


36 ± 2 


Mgl 4571. U9b 


16 


.38 


± 


2 


.98 


136. 


.07 


±3. 


.55 


23 ± 2 


T T n aoc-\ OOT 

tip 4obl.oz7 


19 


.97 


± 





.32 


160. 


.61 


± 0. 


.40 


675 ± 4 


tti„ T , i n o i on o 

le I 4yz4.z9o 


19 


.42 


± 


1 


.46 


130. 


.66 


±2. 


.18 


65 ± 2 


tti„ T , i n o n £?oc 

le I 4yo9.bob 


18 


.22 


± 


2 


.34 


132 


.61 


± 3. 


.18 


29 ± 2 


lei 4957.597 


19 


.55 


± 


1 


.07 


136. 


.19 


± 1. 


.63 


48 ± 1 


tti„ T r nnr iin 

lei 500b. 119 


21 


.48 


± 


2 


.64 


129. 


.10 


± 3. 


.95 


32 ± 2 


le I oUlz.Ubo 


15 


.55 


± 


1 


.77 


140. 


.84 


±2. 


.42 


46 ± 2 


T7i„ T ft n /I 1 < il'7 

be I 5041.447 


20 


.98 


± 


1 


.45 


130. 


.69 


± 1. 


.91 


91 ± 2 


le I 50ol.Do4 


18 


.19 


± 


1 


.64 


139 


.04 


±2. 


.28 


41 ± 2 


T7i T frnTn 

lei 5079.740 


23 


.25 


± 


2 


.15 


129. 


.63 


±2. 


.89 


48 ± 2 


rel oOoo.ooo 


19 


. 10 


± 


1 


::i6 


129. 


.62 


±2. 


.01 


30 i 2 


le I 5107.447 


19 


.03 


± 


1 


.19 


131. 


.36 


± 1. 


.63 


82 ± 2 


TjV. t CT1H /II Q 

re I Oll0.41o 


17 


.30 


± 


1 


.29 


135. 


.01 


± 1. 


.71 


87 i 2 


TV If T C1 OOO 

Mgl 51b7.ozz 


17 


.26 


± 





.45 


139. 


.73 


±0. 


.59 


150 ± 2 


be I 51bo.8y8 


19 


.56 


± 





.75 


131. 


.56 


±0. 


.96 


78 i 2 


TV iT T c 1 TO £J O ,1 

Mgl 517z.bo4 


20 


.00 


± 





.56 


140. 


.22 


±0. 


.78 


254 ± 3 


Mgl 51oo.b04 


17 


.46 


± 





.46 


139. 


.37 


±0. 


.62 


209 i 2 


T71 T coot ion 

lei 52z7.1o9 


21 


.54 


± 





.85 


133. 


.72 


± 1. 


.16 


76 ± 2 


1c I ozby.oo 7 


22 


.13 


± 





.61 


131. 


.14 


±0. 


.85 


173 i 2 


TT1„ T CfOOO noo 

lei 5328.038 


21 


.83 


± 





.12 


132 


.58 


±0. 


.58 


160 ± 2 


T?„ T PT O /I 1 O O ,1 

le I 5o41.zo4 


22 


.18 


± 


1 


.35 


137. 


.92 


± 1. 


.86 


83 ± 2 


tti„ T croTi Hon 

le I 5o71.4o9 


15 


.66 


± 





.60 


130. 


.31 


±0. 


.82 


101 ± 2 


tti„ T cronT "ioo 

lei 5o97.1z8 


17 


.84 


± 





.63 


131. 


.63 


± 0. 


.87 


103 ± 2 


t~^„ T rr a n cr tt /i 

le I 5405.774 


20 


.25 


± 





.85 


135 


.19 


± 1. 


.16 


75 ± 2 


T71 T c ,i on / ■ n rr 

tel 5429.095 


20 


.10 


± 





.60 


133 


.37 


± 0. 


.81 


97 ± 2 


TT1„ T t; A A f OT1 

le I 544b. 871 


19 


.61 


± 





.56 


131. 


.62 


±0. 


. 74 


152 i 2 


T71 T cT/TPrcr /?nn 

le I 5455. b09 


20 


.32 


± 





.56 


130. 


.99 


± 0. 


.76 


121 i 2 


T71„ T CT/inT 

le I 5497.51b 


21 


.74 


± 


1 


.09 


136. 


.31 


± 1. 


.49 


85 ± 2 


le I 5501. 4b5 


18 


.64 


± 


1 


.78 


134. 


.74 


±2. 


.46 


52 ± 3 


T71 T p- trn/^ TTn 

tel 55U0.779 


18 


.37 


± 


1 


.57 


137 


.41 


±2. 


.02 


55 ± 3 


Nai 5889.950 


20 


.63 


± 





.84 


144. 


.01 


± 1. 


.13 


337 ± 7 


Nai 5895.924 


20 


.69 


± 


1 


.02 


144. 


35 


± 1. 


.12 


204 ±7 


Fei 6136.994 


21 


.08 


± 


3 


.63 


133. 


.25 


±4. 


.85 


149 ±7 


Cai 6162.170 


19 


.11 


± 


2 


.74 


140. 


.27 


±4. 


.07 


65 ±6 


Fei 6191.588 


20 


.43 


± 


2 


.43 


139 


.27 


±3. 


.30 


87 ±7 


Fei 6230.723 


20 


.26 


± 


2 


.03 


138. 


.54 


±2. 


.75 


92 ±6 


Fei 6252.555 


16 


.48 


± 


2 


.40 


140. 


.87 


±3. 


.22 


52 ±5 


Fei 6393.601 


21 


.92 


± 


2 


.25 


139 


.69 


±3. 


.10 


65 ±6 


Fei 6400.001 


21 


.00 


± 


2 


.89 


127 


.53 


±3. 


.83 


89 ± 7 


Fe 1 6430.846 


17 


.43 


± 


1 


.97 


139 


.20 


±2. 


.53 


84 ±6 


Fei 6494.980 


18 


.62 


± 


1 


.11 


130 


.75 


± 1. 


.55 


106 ±5 
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Ha 6562.760 18.91 ±0.16 161.66 ±0.20 2894 ± 9 

Fe I 6677.987 23.21 ±2.15 130.50 ±2.70 71 ± 5 

Can 8498.020 20.35±0.13 144.80±0.18 887 ± 4 

Can 8542.090 20.13 ±0.14 144.19 ±0.20 847 ± 4 

Can 8662.140 19.87 ±0.17 142.89 ±0.23 671 ± 4 
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Table A2. Secondary star emission lines in GK Vir. P is the hydrogen Paschen series. After measuring the radial velocity of a line, the 
spectra were shifted to remove the motion and averaged. The Equivalent widths were measured from this averaged spectrum. 



Line 


■7sec 

(km/s) 




^mcas 

(km/s) 




Equivalent 
Width (mA) 


tti /? o to o o rr o 
H16 3703.853 


-45. 


55 


±2. 


07 


200. 


.01 


± 


2 


.93 


68 ± 3 


TTI CT OT1 1 nT1 

H15 3711.971 


-48. 


.49 


± 1. 


83 


203. 


.33 


± 


2 


.49 


81 ± 2 


T oti o o o rr 

lei 3719.935 


-50. 


25 


± 1. 


55 


197. 


.94 


± 


2 


.10 


29 ± 2 


TT1 ^ O TO 1 n /I o 

H14 o721.948 


-46. 


01 


± 1. 


21 


204. 


.48 


± 


1 


.71 


107 ± 2 


TTI Q Q'TQi OTO 

Hid 3734.372 


-49. 


.96 


± 1. 


01 


199 


.31 


± 


1 


.42 


183 ± 2 


TT„ t Of7T7 1 O 1 

lei 3737.131 


-46. 


.50 


± 1. 


00 


199 


.97 


± 


1 


.40 


53 ± 2 


tti o oTrr o 1 co 

H12 o75U.152 


-46. 


68 


±0. 


73 


203. 


.38 


± 


1 


.04 


197 ± 2 


TTI 1 O TTO O A 

Mil 3770.634 


-47. 


06 


± 0. 


.49 


203 


.61 


± 





.70 


239 ± 2 


TTI O O TOT ni o 

H10 3797.910 


-46. 


.32 


±0. 


39 


202 


.72 


± 





.55 


311 ± 2 


T T T OOT O l 7C 1 

He I ool9.7bl 


-48. 


22 


± 0. 


73 


201 


.68 


± 


1 


.05 


69 ± 2 


TTO ooor OOT 

119 3835.397 


-48. 


00 


±0. 


32 


204 


.37 


± 





.46 


331 ± 2 


T nor/? o T 1 

lei 3856.371 


-46. 


.33 


± 1. 


11 


204 


. 18 


± 


1 


.52 


41 i 2 


tto oooo o rr fr 

H8 3889.055 


-47. 


.28 


±0. 


21 


203. 


.86 


± 





.30 


457 ± 2 


o : T oooer rroo 

Di I 3905.523 


-45. 


.96 


± 0. 


69 


205 


.38 


± 





.97 


58 ± 2 


Gall 3933.063 


-46. 


67 


±0. 


28 


205 


.82 


± 





. 39 


208 ± 2 


T T _ on/; TOT 

He I o9b4.727 


-48. 


70 


± 0. 


93 


198. 


.93 


± 


1 


.36 


44 ± 2 


/"""l _ TT onto i r> r\ 

Call oybo.4by 


-46. 


78 


±0. 


31 


206. 


.90 


± 





.40 


158 ± 2 


TT , QO*7n 0.7/1 

He oy(U.U/4 


-47. 


53 


±0. 


20 


205. 


.75 


± 





.27 


529 \ 2 


TT.. _ ^no/ , "1 oo 

Hel 4UZD.loy 


-45. 


77 


±0. 


.49 


204. 


.70 


± 





.70 


77 ± 1 


TTf /I 1 O "1 TO e 

Ho 41Ul,7o5 


-47. 


42 


± 0. 


18 


205. 


.81 


± 





.25 


659 ± 3 


TT i /II /I O TrTO 

He I 414o.759 


-49. 


44 


±0. 


87 


198. 


.33 


± 


1 


.18 


55 ± 2 


/~1 T (no/ 1 TOO 

Cai 422b. 728 


-45. 


78 


± 1. 


44 


209. 


.03 


± 


1 


.95 


36 ± 2 


TT„ T /lOOO TOC 

re I 42o2.72b 


-47. 


77 


± 1. 


15 


201. 


.03 


± 


1 


.55 


41 ± 2 


l^e I 42dd.9d4 


-46. 


93 


± 0. 


92 


198. 


.87 


± 


1 


.26 


54 ± 2 


TT- A o A o A c rr 

H7 4o4U.4do 


-47. 


.48 


±0. 


16 


207. 


.11 


± 





.22 


842 ± 3 


tt„ T ,1 rr 1 z a a 

be 1 4oo 1.544 


-48. 


56 


± 1. 


04 


201. 


.36 


± 


1 


.43 


63 ± 2 


TT-. T A OOT OOO 

He I 4oo7.9z8 


-46. 


72 


±0. 


61 


197. 


.68 


± 





.84 


84 ± 2 


TT„ , /1/1T1 A Of\ 

He I 4471. 4ol) 


-46. 


99 


± 0. 


32 


202. 


.90 


± 





.46 


144 ± 2 


TV T TT /I <01 IOO 

JVLg II 44q1*1oU 


-46. 


51 


±0. 


72 


201 


.59 


± 


1 


.04 


71 ± 2 


he I 4549.467 


-46. 


01 


± 1. 


05 


201 


.43 


± 


1 


.43 


49 ± 2 


Hel 471o.l4b 


-47. 


44 


±0. 


.49 


204. 


.45 


± 





.69 


75 ± 2 


TTO /I Oi? 1 OOT 

Hp 4obl.o27 


-47. 


39 


±0. 


16 


208. 


.67 


± 





.23 


1406 i 4 


tj. T /inoi fion 
He I 4921.929 


-48. 


30 


±0. 


30 


204. 


.02 


± 





.43 


149 i 2 


TT 1 .. TT A (\ O O OOI 

le II 492o.921 


-47. 


35 


±0. 


52 


204. 


.94 


± 





.73 


80 ± 2 


TJ>_ T /IO.tr.7 trO.7 

re I 495 / .59/ 


-46. 


94 


±0. 


98 


205. 


.54 


± 


1 


.41 


59 i 2 


T T T cr o 1 rr i" Trr 

Hel 5U15.o75 


-48. 


93 


±0. 


26 


205. 


.02 


± 





.37 


197 ± 2 


TT^ tt Cfl1 O AO A 

r C II 5Ulo.4o4 


-48. 


60 


±0. 


50 


203. 


.93 


± 





.73 


105 i 2 


HiT T FT 1 T OOO 

JVLg I 51b7.o22 


-49. 


24 


±0. 


63 


205. 


.19 


± 





.90 


90 ± 2 


re I 51bo.o9o 


-45. 


00 


±0. 


51 


207. 


.13 


± 





.73 


112 i 2 


TV T T rr 1 TO CO/I 

JVLg I 5172. bo4 


-46. 


67 


±0. 


96 


207. 


.01 


± 


1 


.41 


127 ± 2 


IVLg I 51oo.uU4 


-48. 


44 


±0. 


66 


205. 


.68 


± 





.99 


136 i 3 


Tn„ T rr oo t 1 00 

lc I 5227.189 


-46. 


44 


±0. 


76 


200. 


.01 


± 


1 


.10 


70 ± 2 


T7< T rro/"o rroT 

re I 5269.537 


-48. 


.88 


± 0. 


85 


198. 


.55 


± 


1 


.27 


92 ± 2 


TT,. TT tTOTj^ OOO 

re II 527b. UU2 


-48. 


15 


±0. 


77 


202. 


.46 


± 


1 


.11 


80 ± 2 


TT 1 .. TT rroi£? £M rr 

le II 5316.615 


-46. 


43 


± 0. 


57 


201. 


.34 


± 





.82 


102 i 2 


T7< T rTOOO OOO 

1c I 5328.038 


-47. 


10 


±0. 


67 


203. 


.37 


± 





.99 


126 i 2 


tt 1 ^. T rr /i rr tt a 

re I 54U5.774 


-48. 


88 


± 1. 


95 


205. 


.80 


± 


1 


.59 


51 ± 2 


Tn„ T rr .1 oo /^orr 

re I 5429. o95 


-46. 


80 


±0. 


96 


207. 


.31 


± 


1 


.43 


73 ± 2 


tt„ T rz a d a rro/i 

re I 54o4.524 


-48. 


12 


± 2. 


77 


198. 


.15 


± 


3 


.37 


25 ± 2 


T7< t- A A C O T1 

le I 544b. 871 


—45. 


26 


± 0. 


97 


206. 


.78 


± 


1 


.41 


53 ± 2 


tt„ T rr a rr rr r*oo 

l^e 1 5455. bU9 


-46. 


39 


±0. 


82 


204. 


.71 


± 


1 


.20 


62 ± 2 


tt.. 1 rroTrr a i o 

He I 5875.618 


-47. 


07 


±0. 


29 


207. 


.82 


± 





.41 


623 ± 7 


Nai 5889.950 


-48. 


81 


±1. 


07 


211. 


.58 


± 


1 


.49 


160 ± 8 


Nai 5895.924 


-49. 


29 


±1. 


13 


209. 


.46 


± 


1 


.59 


124 ±8 


Fei 6136.994 


-45. 


82 


±3. 


19 


207. 


.64 


± 


3 


.19 


112 ±8 


Cai 6162.170 


-49. 


04 


±2. 


86 


210. 


.47 


± 


2 


.27 


31 ± 5 


Fei 6191.588 


-48. 


29 


±1. 


77 


205. 


.30 


± 


2 


.30 


47 ±7 


Mgn 6346.962 


-45. 


17 


±0. 


89 


203. 


.44 


± 


1 


.21 


123 ±6 


Fei 6393.601 


-46. 


.88 


±2. 


70 


198. 


.16 


± 


3 


.22 


50 ±7 


Fei 6400.001 


-46. 


92 


±1. 


74 


208. 


.27 


± 


2 


.57 


57 ±6 


Fen 6456.383 


-46. 


43 


±1. 


70 


204. 


.89 


± 


2 


.21 


125 ± 7 


Cai 6462.567 


-47. 


03 


±1. 


90 


209. 


.01 


± 


2 


.58 


37 ±5 


Fei 6494.980 


-46. 


55 


±1. 


30 


209. 


.76 


± 


1 


.74 


56 ± 5 
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Precise parameters for SDSS J1212-0123 and GK Vir 21 



Ha 6562.760 


-46 


.77 


± 


0. 


.25 


210 


.73 


± 





.33 


4073 ± 13 


He I 6678.149 


-47. 


.50 


± 





.51 


206 


.26 


± 





.61 


645 ± 7 


He I 7065.188 


-46. 


.25 


± 





.54 


205 


.69 


± 





.64 


820 ± 7 


He I 7281.349 


-48. 


.65 


± 


0. 


.64 


202 


.81 


± 





.82 


344 ± 8 


Oi 7771.944 


-47. 


.61 


± 


0. 


.61 


203 


.64 


± 





.81 


223 ± 5 


Oi 7774.166 


-48. 


.06 


± 


1. 


.39 


202 


.81 


± 


1. 


.12 


234 ± 5 


Oi 7775.388 


-46. 


.94 


± 


1. 


.55 


204 


.89 


± 


1. 


.38 


131 ± 5 


Mgn 7896.368 


-48. 


.23 


± 





.90 


202 


.55 


± 


1. 


.30 


131 ± 7 


Nai 8183.256 


-49. 


.93 


± 


1 


.08 


206 


.13 


± 


1. 


.62 


74 ± 6 


Nai 8194.824 


-46. 


.10 


± 


1. 


. 16 


206 


.71 


± 


2 


.13 


89 ± 7 


Fei 8327.056 


-47. 


.02 


± 


1. 


.04 


200 


.82 


± 


1. 


.54 


158 ± 8 


Ci 8335.150 


-46. 


.96 


± 





.69 


207 


.21 


± 





.94 


283 ± 7 


O I 8446.359 


-47. 


.51 


± 


0. 


.73 


206 


.06 


± 





.98 


340 ± 9 


Can 8498.020 


-46. 


.89 


± 





.19 


207 


.18 


± 


0. 


.27 


925 ± 9 


Can 8542.090 


-47. 


.39 


± 





.15 


207 


.64 


± 





.21 


1322 ± 9 


P14 8598.392 


-48. 


.88 


± 


2. 


.00 


200 


.56 


± 


2 


.78 


256 ± 12 


Can 8662.140 


-47. 


.65 


± 


0. 


.19 


207 


.95 


± 





.29 


901 ± 9 


P13 8665.019 


-45. 


.51 


± 


1. 


.88 


202 


.98 


± 


2 


.68 


489 ± 12 


Fei 8688.625 


-48. 


.01 


± 


1. 


.07 


203 


.53 


± 


1. 


.66 


176 ± 9 


P12 8750.473 


-44. 


.68 


± 


1. 


.71 


203 


.42 


± 


2. 


.50 


719 ± 13 


Mgi 8806.757 


-46. 


.67 


± 


0. 


.72 


205 


.38 


± 





.99 


244 ± 8 


Pll 8862.784 


-45. 


.69 


± 


1. 


.35 


202 


.23 


± 


1. 


.93 


1010 ± 13 


Can 8912.070 


-49. 


.99 


± 





.78 


202 


.51 


± 


1. 


.19 


194 ± 9 


Can 8927.360 


-49. 


.91 


± 


0. 


.76 


202 


.22 


± 


1. 


.10 


207 ± 9 


P10 9014.911 


-47. 


.08 


± 


0. 


.66 


201 


.91 


± 





.95 


1337 ± 13 


Ci 9061.430 


-46. 


.69 


± 





.88 


202 


.20 


± 


1. 


.12 


638 ± 10 


Ci 9078.280 


-46. 


.30 


± 





.85 


203 


.85 


± 


1. 


.08 


255 ± 9 


Ci 9088.510 


-47. 


.30 


± 





.79 


208 


.63 


± 


1. 


.07 


378 ± 10 


C I 9094.830 


-45. 


.99 


± 


0. 


.90 


201 


.89 


± 


1. 


.19 


385 ± 10 


Ci 9111.800 


-48. 


.71 


± 


0. 


.79 


206 


.42 


± 


1. 


.01 


281 ± 9 


Fei 9212.981 


-48. 


.42 


± 


1. 


.14 


207 


.64 


± 


1. 


.67 


212 ± 12 


Mgn 9218.248 


—46. 


.67 


± 





.48 


204 


.95 


± 





.70 


493 ± 12 


P9 9229.015 


-49. 


.80 


± 


0. 


.54 


203 


.76 


± 





.78 


2079 ± 20 


Mg n 9244.266 


-46. 


.71 


± 


0. 


.63 


203 


.53 


± 





.91 


289 ± 11 


Ci 9405.730 


-47. 


.00 


± 


0. 


.77 


208 


.83 


± 


1. 


.01 


799 ± 19 


Pe 9545.972 


-47. 


.40 


± 





.43 


205 


.84 


± 





.62 


2825 ± 26 


Ci 9620.800 


-48. 


.95 


± 


1. 


.00 


204 


.99 


± 


1. 


.11 


243 ± 15 


Ci 9658.440 


-45. 


.49 


± 


0. 


.99 


208 


.87 


± 


1. 


.10 


202 ± 18 


P5 10049.37 


-48. 


.12 


± 





.72 


206 


.70 


± 





.65 


3576 ± 60 
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